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ABSTRACT 

We perform a detail analysis of two models of neutrino CC A production on free nucleons. First model 
is a standard one based on nucleon-Delta transition current with several form-factors. Second model is a 
starting point for a construction of Marteau model with sophisticated analytical computations of nuclear 
effects. We conclude that both models lead to similar results. 
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SHORT TITLE: On two weak CC A production models 



1. INTRODUCTION 

Future precise neutrino measurements e.g. of 
f9 13 require better understanding of neutrino in- 
teractions in few GeV region pQ. In general, in 
description of neutrino-nucleon interaction ver- 
tex three processes are distinguished: quasielas- 
tic, single pion production (resonance excitation 
region) and more inelastic processes taken into 
account by the formalism of deep inelastic scat- 
tering. Weak single pion production is therefore 
a part of more complicated dynamics. It gives an 
important contribution to cross section in 1 GeV 
region and has to be treated with care. In the past 
it was a subject of many theoretical studies A 
sample of existing experimental data is not con- 
clusive as measurements were made with a typical 
precision of 20-25% [Hj- From a point of view of 
Monte Carlo simulation codes there seems to be 
an agreement that Rein-Sehgal 0J model is most 
reliable. It includes contributions from 18 reso- 
nances with masses up to 2 GeV, their interfer- 
ence terms together with a non-resonance back- 
ground. Recent developments in quark-hadron 
duality suggest however that there is no need to 
consider so many resonances: contributions from 
most of them can be described in average by suit- 
ably modified PDF's [5]. When reaction takes 
place on nuclear targets resonance contributions 
are additionally smeared out by Fermi motion. A 
conclusion is that probably only the A excitation 
has to be treated independently 0. 

One way to describe A excitation is to con- 
struct a current < AlJ^I^V > with phenomeno- 
logical form-factor constrained by CVC and 
PC AC arguments 0. There have been also at- 
tempts to calculate such form-factors from first 
principles in the quark model 8 . A precision 
with which form-factors are known cannot be bet- 
ter then experimental uncertainties. 

Few authors tried to discuss nuclear effects in 
single pion production in a framework of more 



systematic theoretical schemes. One of such mod- 
els was developed by Marteau [Hj. It includes: 
Fermi motion, Pauli blocking, elementary lp-lh, 
lA-lh and 2p-2h excitations, modification of A 
width in a nuclear matter, RPA corrections and 
finite volume effects. Before all the nuclear ef- 
fects are taken into account, the model is based 
on simplified dynamical assumptions about CC 
neutrino A excitation on free nucleons. These 
simplifications are necessary in order to perforin 
calculations of nuclear effects in compact and el- 
egant way. 
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Figure 1. Total cross section for ir + produc- 
tion via A ++ excitation. Experimental points are 
taken from |3j . No constraints on hadronic invari- 
ant mass are imposed. 



The Marteau model has been discussed in both 
experimental and theoretical context. On exper- 
imental side JOI the model is used in the Monte 
Carlo code of K2K collaboration to give predic- 
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Figure 2. Total cross section of tt + production 
via A ++ excitation with a restriction on invari- 
ant hadronic mass W < lAGeV . Experimental 
points are taken from [2j . 



Figure 3. Total cross section of ir + production 
via A ++ excitation with a restriction on invari- 
ant hadronic mass W < l.QGeV. Experimental 
points are taken from . 



tions for NC 7r° production. There is an open 
dispute as its predictions are much smaller then 
those of Rcin-Scghal model for coherent pion pro- 
duction . It has been also investigated for the 
problem of pion-less A decays in the nuclear mat- 
ter UJ. The Marteau model was mentioned by 
several authors dealing with theoretical descrip- 
tion of nuclear effects in neutrino interactions in 
GeV region [H] . 

Due to the lack of experimental data it is diffi- 
cult to evaluate how successful the model is in de- 
scribing nuclear effects. Certainly it is important 
to investigate its properties by performing com- 
parisons with other approaches. An important 
test is to investigate basic dynamical assumptions 
of the model and compare results with predictions 
of other models and with existing experimental 
data on free nucleons. Such tests are presented 
in this paper. Our main conclusion is that CC 
neutrino-/ree nucleon interaction model being a 
basis for the Marteau model (we call that model 
pre-Marteau model) leads to very similar results 
as the model |7]. 



2. RESULTS 

Marteau model [H] is a sophisticated model con- 
taining in one theoretical scheme neutrino CC 
quasielastic and A excitation reactions on nuclei. 
We present below its basic assumption for the dif- 
ferential cross section for CC neutrino A ++ exci- 
tation on free nucleoli: 



d 2 <r 



A++ 

M 



G 2 cos 2 e c qMlT A 
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In order to avoid a confusion we call this model 
pre-Marteau model. The above formula is derived 
in Appendix A, also the notation is explained 
there in detail. 

In this paper we compare predictions of pre- 
Marteau model with the model for neutrino CC 
A excitation based on nucleon-A transition cur- 
rent with several form-factors. In the latter model 
[7] (we call it Form-Factor Model) one calculates 
cross section in a standard way. Straightforward 
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Figure 4. Total cross section of 7r + production via 
A ++ excitation with a decomposition into longi- 
tudinal and transverse contributions. 



Figure 5. Differential cross section of invariant 
hadronic mass for ir + production via A ++ exci- 
tation. Neutrino energy is E v = lGeV. 



computations lead to expression 
a <7 FF 



G 2 cos 2 9 C T A L^Wqdwdq 
6Att 2 E 2 (W -Ml) +r^/4 



(2) 



In order to calculate Ti.^" one introduces 
nucleon-A transition current [15| 



M 



(g^q - q»r) 



M 2 

c£{Q 2 ) 



Vvp'-«y a ))75+ 



M 2 



(g^q-p'-q^p' a ) 



(3) 



+C£(Q 2 )g^ + ^f±q»q a )u(p). 

where q^ = p^ — p'^. A factor yo is present in 
the current for A ++ production. In our numerical 
calculations we used the form factors from |13j . 
Equations (1) and (2) as they are written very 
similar. In fact the kinematics is the same in both 
cases. The inequivalence of two expressions comes 
from hadronic tensors (i? A " v )^ =0 1V 1V which are 
calculated in different ways. 



A comparison of predictions of two models was 
done for a free nucleon target because the existing 
experimental data applies to this situation. We 
compared total cross section in the energy range 
up to 5 GeV with data from [j^. We presented 
plots without bound on the invariant hadronic 
mass and with bounds 1.4 GeV and 1.6 GeV as 
such experimental data is available. 

We made a decomposition of the total cross sec- 
tion into two parts: "transverse" and "longitudi- 
nal" according to spin-isospin operators present 
in the hadronic current. In the pre-Marteau 
model longitudinal operators are present only in 
Hqq, H33, i?03 and transverse only in -ffn, Hyi (in 
the frame q = (0, 0, q)). Therefore we define "lon- 
gitudinal" and "transverse" parts as those coming 
from corresponding terms in L^ V H^ V . 

For two theoretical models we performed also a 
comparison of differential cross sections for neu- 
trino energy of IGeV and 2GeV of: invariant 
mass, energy transfer and cos (5a, an angle be- 
tween momenta of incident neutrino and A. 

Our conclusions are the following: 
i) For neutrino energy lower than E v — 5GeV 
(Fig. 1) pre-Marteau model gives rise to higher 



5 




Figure 6. Differential cross section of energy 
transfer for n + production via A ++ excitation. 
Neutrino energy is E v = lGeV. 



Figure 7. Differential cross section of energy 
transfer for n + production via A ++ excitation. 
Neutrino energy is E v = 2GeV. 



values of the total cross section but both models 
approximately agree with the existing experimen- 
tal data. 

ii) For E v = 5GeV Form- Factor Model predicts 
the values of the total cross section greater than 
pre-Marteau model. Pre-Marteau model cross 
section becomes approximately constant while 
Form-Factor Model cross section still increases 
with neutrino energy. 

iii) With constraints on the value of invariant 
hadronic mass (Figs 2, 3) both models give rise 
to predictions similar in shape and in agreement 
with experimental data. Values of total cross sec- 
tion of pre-Marteau model are about 20% higher 
than predicted by Form-factor model. 

iv) Contributions from "longitudinal" and 
"transverse" contributions (Fig. 4) to the total 
cross section are similar in both models. 

v) Invariant hadronic mass distributions (Fig. 
5, 6) are very similar, the only difference is in 
scale and comes from different values of the total 
cross section. For higher values of neutrino energy 
Form- Factor Model gives rise to greater contribu- 
tion from higher values of W. This is why cuts 
on the invariant mass are more restrictive for that 



model. 

vi) Differential cross sections of energy transfer 
(Fig. 7, 8) for both models are similar in shape. 
There is an surprising decline of differential cross 
section in energy transfer at about v = Q.7GeV 
for E v = IGeV (Fig. 7) and at about v = l.7GeV 
for E v — 2GeV (Fig. 8) present in predictions of 
both models. A possible explanation is kinemat- 
ical in origin. For a fixed value of energy transfer 
v the integration domain in momentum transfer 
q is 

qe I y/v 2 -m 2 + 2E(E -v)-p, 



min(y^ 2 - m 2 + 2E(E - v) + p, ( 4 ) 
y/v 2 -ml + 2M (iz-m^^j. 

where 

p = 2E x /(E-v) 2 -m 2 . 

It turns out that at above mentioned values of 
v and E the allowed kinematical domain begins 
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Figure 8. Differential cross section of cos 9 a for 
7r + production via A ++ excitation. Neutrino en- 
ergy is E v = lGeV. 



Figure 9. Differential cross section of cos 9 a for 
7r + production via A ++ excitation. Neutrino en- 
ergy is E v = 2GeV. 



to decrease quickly with v because the first ar- 
gument in min becomes smaller then the second 
one. 

vii) In the differential cross section of cos 0a 
(0 A is angle between A and incident neutrino mo- 
menta) (Fig. 9, 10) one can notice maxima at 
values cos (5 A ~ 0.875 (0 A ~ 29°) for E v = IGeV 
(Fig. 9) and cos0 A ~ 0.775 (0 A ~ 39°) for 
E v = 2GeV (Fig. 10). The maxima are present in 
predictions of both models. For neutrino energy 
E v = 2GeV (Fig. 12) there is is also a second 
maximum in the forward direction. The shape 
of differential cross-sections is similar to the one 
derived in [T5| . 

Our final conclusion is that pre-Marteau model 
for A ++ excitation leads to close to standard be- 
havior of 7r + production cross section and it is le- 
gitimate to use it in sophisticated computations 
of nuclear effects and in MC codes . 



APPENDIX A 

A logic of the Marteau model can be under- 
stood if one starts from the differential cross sec- 



tion for quasi-elastic process n — > /J, p in the 
Fermi gas model |14j : 

d 2 af l G = J d 3 pdvdq8{k F - |pl)0(|p + q\ - k F ) 

xSiE + Ef-Ep-E?) (5) 
G 2 cos 2 6 c N kF qMM' 



l&EpEp-KE 2 



-Lni/H^". 



fc M , p^, p ,fl denote 4-momenta of: neutrino, 
charged lepton, target and recoil nucleons. M 
and M' are masses of target and ejected nucleons. 
In the case of quasielastic process they are taken 
as equal, q^ = fc' 1 — k'^ = {u, q) is energy and 
momentum transfer. and H^ v are leptonic 

and hadronic tensors: 



(6) 



If 11 " = -r-^TrfP'tf + M)t»(p' + M')J , (7) 

(8) 



1 
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f v = 7 °(r' i )t 7 , 
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T^ = F 1 {Q 2 ) 1 ^iF 2 {Q 2 )a^^ 
+G A (Q 2 )n 5 + G P (Q 2 ) l5 £ I . 



(9) 



Fx, F2, Ga and Gp are standard form- factors 
HS| N kF = g|pr = j^g. We assume nucleus 
of atomic number A to contain equal numbers of 
protons and neutrons. 

In a good approximation (since we will actually 
calculate cross section for a reaction on free nu- 
cleons our derivation becomes exact; we perform 
all the steps of computations starting from the 
Fermi gas in order to check normalization fac- 
tors) in E } E - we put p — (thus p' = q). 
The dependence on p factorizes and we define 



/ 



d 3 P 6(k F -\p\)9(\p + q\-k F ) 
x5(E + Ep- Ep ~Ep 
V 



(10) 



2N kF TT 



TmIL~N-h(v,q)- 



2 qel _ 
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FG - -2mYl N ^ h (v,q) 
G 2 cos 2 6 c M'qV 



(11) 



32ir 2 E 2 E ir 



L^(H^)p =0 dudq 



In the original Marteau approach non- 
relativistic nucleon's kinematics is used and 
Tm%_/,(!/,g) is a Lindhard function, the 
particle-hole polarization tensor, an object which 
accounts for Fermi motion and Pauli blocking |17| 
(in [§] higher order corrections in are consid- 
ered). 

In the limit p - 
position of 



A I 



using non-relativistic decom- 



terms coming from different spin operators. For 
example we calculated |16| : 



J° = Njtf, ( l(Fx - F 2 ±- ) 

pVs \ v 1 2 M' + E p , 



-a ■ q- 



JG A - vG P 
M' + E p , 



(14) 



where N p < — y Ev '^f and </>'s describe non- 
relativistic spinors. 

Marteau A excitation model is defined by |T5) : 



(i) substitution M 1 = M A ; 

(ii) multiplication of form-factors by the nu- 
merical factor 4.78 = (#^) 2 ; 

(iii) elimination of " charge" terms (spin opera- 
tor is necessary to produce a particle of spin 3/2) 
from H^ v - we call the new tensor H^ v with the 
numerical factor 4.78 included in its definition; 

(iv) substitution of Ipy-h by Ha— h the polar- 
ization tensor for A -hole excitation: 



The expression for the cross section takes form ImIL A -h{v, q) 



32 1 

m 2 a t a 



d 3 P 0(k F - \p\) 



(15) 



{W 2 - M A ) 2 + M A Tf 



(the factor comes from summation over 



isospin and spin degrees of freedom), 
W 2 = {Ep + vf-{p + qf- 



(16) 



(v) inclusion of RPA correlations and local den- 
sity effects. 

In our derivation we restricted ourselves to 
steps (i-iv) and we obtained a model for A ex- 
citation on nuclei in the free Fermi gas approxi- 
mation: 



J" = u{q)r>*u{d). 
one can identify in 



(12) 



(13) 



A 1 ^ 



G 2 cos 2 6 c VM A q 
32-K 2 E 2 Ef 



(17) 



XImIlA-/,(y, q)L )LV {H ilv )p^dvdq. 



In the limit kp — > we obtained pre-Marteau 
model (in this limit target nucleon is at rest and 
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our evaluation of 



becomes exact) 
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Af 2 r A 



(18) 



((Af + - g 2 - Afi)' + M^r^ 
Finally (after fcp's get properly cancelled) 
A _ G 2 cos 2 6> c Af A r A A? 



d 2 CT 



36tt 2 E 2 (M + v) 



(19) 



((M + //) 2 -q 2 



Mi) 



M 2 ri 



r A is defined as (r = 115MeV) 

lA=l0 7T7 \o "TT7i (20) 

where q cm (W) is the pion momentum in A (of 
mass W) rest frame. 

Pre-Marteau model provides a prediction for an 
overall A production i.e. for a sum over isospin 
degree of freedom. Without nuclear effects rel- 
ative probabilities to produce isospin states is 
given by a ratio of Clebsh-Gordan coefficients. 
Thus for neutrino induced reaction the probabil- 
ity to produce A ++ is three times the probability 
to produce A + . In this paper we present a com- 
parison for A ++ production. It is because in the 
measurements of the invariant hadronic mass dis- 
tribution for the process v^n — > fi~pTr + there is 
a sharp resonance peak at W ~ 1.2 GeV while 
in the channels v^p — > /i~?i7r + and v^p — > n~pir° 
peaks are smeared out It is clear that cor- 
rect description of the last two channels requires 
an addition of non- resonant contribution or/and 
contributions from other resonances while in the 
first one A ++ production cross section can be 
meaningfully compared with the data. Prediction 
for A ++ production per nucleon is thus obtained 
by dividing by 4 anc ^ multiplying by | . The 
final formula for A ++ excitation cross section per 
proton in the pre-Marteau model reads: 



M 



_,2_ A ++ _ G 2 cos2 O c qM\T A 
24tt 2 E 2 (M + v) 

L^{H^)j; =0 dvdq 

((Af + v) 2 - q 2 - M 2 A ) 2 + M|r A 



(21) 
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